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We study the interplay between charge doping and intermolecular distance in the polymerization 
of C60 fullerene chains by means of density functional theory (DFT)-based first principle calculations. 
The potential energy surface analysis shows that both the equilibrium intermolecular distance of 
the unpolymerized system and the polymerization energy barrier are inversely proportional to the 
electron doping of the system. We analyze the origin of this charge-induced polymerization effect 
by studying the behavior of the system's wavefunctions around the Fermi level and the structural 
modifications of the molecules as a function of two variables: the distance between the centers of 
the molecules and the number of electrons added to the system. 
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I. INTRODUCTION 



The formation of covalent intermolecular bonds in solid 
Ceo has been observed to occur under a great variety of 
chemical and physical conditions. Experimentally, the 
addition reaction between neighboring fullerenes was first 
observed upon visible or ultraviolet light treatment [l[ to 
crystals of pristine Cgo- Later, the application of high 
pressure and high temperature conditions as well 

as electric field stimulated techniques [f|, have been re- 
ported as efficient methods for the formation of chains 
of polymerized fullerene molecules in solid Cgo- The in- 
tercalation of alkali metal atoms in fullerite crystals has 
been observed to favor or prevent the formation of Ceo 
intermolecular covalent bonds, depending on the nature 
of the intercalated atoms and on the stoichiometry of 
the crystal. Low content or low size of alkali atoms ac- 
tually favor the formation of Cgo polymers. In particu- 
lar, the MCeo compounds (M=K, Rb and Cs) undergo a 
structural transition from a high-temperature monomer 
to a polymer phase for T<350 K [Q-Q. The Cgo-based 
compounds intercalated with light alkali atoms (Li and 
Na), like Na2RbCeo, Na4C6o and Li4Cgo have been ob- 
served to spontaneously form ID and 2D polymers, re- 
spectively, at ambient conditions [ToL IT]] ] . The polymeri- 
azation between two neighboring molecules usually takes 
place via a [2+2]-cycloaddition. In fact, the forbiddncss 
of the suprafacial [2+2]-cycloaddition, as stated by the 
Woodward-Hoffmann rule, is removed by the symme- 
try breaking of the orbital spectrum resulting from elec- 
tron excitation, electron doping or a structural displace- 
ment occurring in the transition state [l2| • On the other 
hand, fullerides with high stoichiometries (e.g. MgCgo, 
M=alkali metal) of light or heavy alkali atoms have not 
been observed to form Ceo polymers at ambient condi- 
tions of temperature and pressure. We recently showed 
that for high content of Rb or Cs metal atoms a signif- 
icant smaller distance between molecules, compared to 



pristine Cgo, is required to observe a transition to a poly- 
merized phase jl3l - il5| . All this suggests that the inter- 
play between electron doping and lattice parameter is es- 
sential to understand the polymerization. Theoretical in- 
vestigation about the stability of negatively charged iso- 
lated Ceo molecules [HI and of anionic fullerenes chains 
connected through face-to-face he xag ons have been al- 
ready reported in earlier works 17]. Also, the elec- 



tronic properties, geometric structure and optical prop- 
erties of negatively charged dimers have been studied 
by employing semiempirical methods [HI, [l!| . More re- 
cently, the dimerization of Cgo has been studied in terms 
of the intermolecular interaction potential between iso- 
lated molecules by taking into account inter-ball donor- 
acceptor contributions [201 ] . In the present work, we per- 
form a detailed study of the interplay between electron 
doping and inter-Cgo distance on the formation of inter- 
ball covalent bonds in attempt to understand the exper- 
imental results and guide future works. This knowledge 
is of importance to understand the interplay of charge 
transfer and steric effects on the polymerization of Cgo 
which can play a fundamental role in fullerene engineer- 
ing of new materials for instance through high pressure 
and high temperature |21| . This work reports a first prin- 



ciples calculations study based on the density functional 
theroy (DFT) on a ID Cgo chain with explict addition 
electrons. 



II. METHODOLODY 

Calculations arc performed with the SIESTA [22[ 
method, where the eigenstates of the Kohn-Sham hamil- 
tonian are expressed as linear combinations of numerical 
atomic orbitals. Many body effects are described within 
the local density approximation (LDA) to the exchange 
and correlation potential [23j | . We have used a varia- 
tionally optimized (2(| H?} double-^ polarized basis sets. 
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Real space integrals were performed on a mesh with a 
250 Ry cutoff. Core electrons are replaced by nonlocal, 
norm-conserving fully separable Trouiller-Martins pseu- 
dopotentials. In the calculations 2s and 2p electrons of C 
atoms were explicitly included in the valence. For each 
calculation, we have minimized the total energy of the 
system as a function of the atomic coordinates imposing 
a tolerance in the largest component of the forces of 0.03 
eV/A. A geometry optimization is performed for each 
intermolecular distance and for each charged state. 

Our study has been performed as a function of two 
variables: (i) the distance R (see Figure [1]) between the 
centers of the molecules along the chain (uniformly var- 
ied from 8.50 A to 10.85 A); (ii) the added charge q, 
consisting on the addition of 1, 2, 3 and 4 electrons to 
the system. A neutralizing jcllium background of oppo- 
site charge ensures system neutrality as a whole in such 
a way that ptot{r) = p(r) + Q, where ptot{r) is the elec- 
tronic density of the system with the extra charge and 
Q is the constant density background neutralizing the 
system [H,[29|]. 

We employed a cubic supercell (of lattice parameter 
equal to L=2R), containing two molecules along one of 
the Cartesian axis in order to allow their relative reori- 
entation during the study. In this way, since the two 
nearest molecules at R distance are placed along the 
<001> direction, the distance between two neighboring 
molecules along the <010> and the <100> directions is 
large enough to allow us to neglect their interaction. We 
therefore refer to the studied systems as fullerene chains. 



mixed sp 2 -sp 3 hybridized bond. 



III. RESULTS AND DISCUSSION 

A. Potential Energy surface 

The total energy of the system as a function of the 
intermolecular distance R and extra charge is reported 
in Figure The total energy for the charged systems 
has been corrected in order to remove the spurious con- 
tribution due to the periodic boundary conditions [3l[ as 
explained below. 





R 



FIG. 1: Arrangement of the molecules in the unit cell. R is 
the center-to-center distance between neighboring fullerenes. 
Atoms numbered as 1, 2, 3 and 4 participate in the [2+2]- 
cycloaddition. 

The starting relative-orientation of the molecules cor- 
responds to that most frequently observed in the one- 
dimensional polymerization of solid fee Cgo [3T.[30|. where 
a double bond of one molecule (atoms 1-2 in Figure[T]) op- 
poses a double bond of the neighboring molecule (atoms 
4-5 in Figure [1]). In this way, the polymerization takes 
place via a [2+2]-cycloaddition mechanism, meaning that 
two double bonds of two adjacent molecules break up in 
order to form two covalcnt bonds between the molecules. 
In the so formed four membered ring, each carbon atom 
is covalently bonded to other four C atoms through a 



FIG. 2: Total energy profile as a function of the intermolecular 
distance for different anionic states. Curves are shifted in 
such a way that the minimum energy value for each curve 
is set to zero. The energy values for points at R > R po i 
have been corrected by the Madelung term. The inset shows 
the decrease with intermolecular distance of the square root 
of the second moment of the charge density distribution, a, 
indicating a more localized charge distribution around the 
fullerene. 

The first and second minima of the total energy profile 
correspond to polymerized (at R=R po {) and unpolymer- 
ized (at larger R) molecules, respectively, as shown in the 
following sections. For all the studied anionic states, as 
well as for the neutral case, the most stable configuration 
is the polymerized one. The polymerization at R < R po i 
is prevented by an energy barrier which is strongly de- 
pendent on the charge of the system. 

The sharp energy jump observed in Figure [3] at dis- 
tance R ~ Rpoi is a clear indication of the avoided cross- 
ing of the two adiabatic surfaces (one for two non bonded 
Ceo and the other for the polymerized molecules). Since 
the transition state of the polymerization is never ob- 
served, the energy barriers here reported are related but 
do not correspond to actual activation energies. 

In the region of R > R po i we can approximate our 
system as a purely molecular system and we apply the 
Madelung energy correction for a lattice of point charges 
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FIG. 3: Electron affinities for the chain of unpolymerized 
molecules as a function of the intermolecular distance R. 
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immersed in a neutralizing jellium [3l|. The energy is 
then expressed by the following equation: 

E(L) = E LDA -^ + 0(L~ 3 ) (1) 

where q is the extra charge and a is the Madelung con- 
stant whose value is known for the simple cubic geom- 
etry. At R < Rp i, although the ground state of the 
system is that of a polymerized phase, this correction is 
still a reasonable approximation since the localization of 
the charge between the Cgo molecules is negligible com- 
pared to the charge around the fullerenes. Also, no dis- 
continuity of the localization of the charge density on the 
fullerene is observed at R=R po i as shown in the inset of 
Figure [2j The measure of localization of the electronic 
charge on the fullerene is made by computing the spher- 
ical second moment of the charge density distribution, 
i.e. f(r— {r)) 2 p(r)dr. We observe that the total energy 
of the system drastically decreases upon addition of ex- 
tra charge indicating an extremely high positive electron 
affinity (EA). In Figure 2] we report the EAs (from the 
first-EA to the fourth-EA) for a linear chain of molecules 
as a function of the intermolecular distance, for R > R po i. 
The (first) electron affinity for the two molecule-system 
is higher than the measured values for the isolated Ceo 
(2.67 eV HH) and similarly computed values (2.75 eV 
obtained with GGA by Pederson et ai. [H|). This is con- 
sistent with the decrease in quantum confinement of the 
extra electron in the two molecule-system as compared 
to the isolated fullerene. The monotonic decrease of the 
electron affinity as a function of the inter-ball distance 
can be associated as well to an increase in quantum con- 
finement. The delocalization of the HOMO wavefunction 
decreases with increasing intermolecular distance (inset 
of Figure [2]). The second momentum of the charge den- 
sity distribution at #=10.85 A (-0.3 A) is 30% smaller 
than at i?=9.05 A (~0.4 A). 

The binding energy between neutral Cgo molecules as 
obtained by the total energy curve of Figure [2] is 0.15 



FIG. 4: Top panel: decrease of the energy barrier for the 
[2+2]-cycloaddition reaction as a function of the charge. Bot- 
tom panel: decrease of the position of the second minimum 
of the total energy curve as a function of the charge. 



eV/2 molecules and the equilibrium distance around 10.2 
A. These values are very similar to those obtained from 
the Lennard- Jones potential (0.26 eV/2 molecules and 
10 A, respectively) reported in Ref. [33| and computed to 
match ex per imental values of cohesive energy and lattice 
constant [34j. Such agreement can be attributed to the 
LDA ovcrestimation of binding which implicitly accounts 
for dispersive forces. 

While the equilibrium interaction distance R po i slightly 
changes as a function of number of electrons added to the 
system, the energy barrier height is significantly affected 
by the different charge state. The existence of such an 
energy barrier for the formation of covalent bonds be- 
tween fullerenes has been already observed experimen- 
tally [HI, HH and theoretically studied [2(| ■ The value of 
the energy barrier obtained for the different anionic states 
is reported in Figure [3J For the neutral system this value 
(2.64 eV) significantly exceeds the activation energy ob- 
served experimentally for photostimulated dimcrization 
(1.25 eV [13) and pressure- induced polymerization in 
pristine Ceo_Q.40 eV [HI) and also the calculated one 
(1.49 eV in |2fJ). This is due to the fact that the energy 
barriers we show only represent an upper limit to the ac- 
tivation energy. We observe that upon electron doping 
the energy barrier dramatically decreases from a maxi- 
mum of 2.64 eV for the neutral system down to a mini- 
mum of 0.67 eV for a system with 4 additional electrons. 
This result proves that the charge actually determines 
the thermodynamic conditions at which polymerization 
takes place. Another important effect of the charge is 
to decrease the equilibrium intermolecular separation for 
the unpolymerized phase, as seen from the charge depen- 
dent position of second minimum of the potential energy 
reported in Figure OH We observe a 10% change in the 
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equilibrium distance between the neutral and q=A case 
due to the progressive population of bonding states with 
charge as discussed in the following section. 

B. Analysis of the wavefunctions close to the Fermi 
energy 

The Cgo molecules in the polymeric chain belong to 
symmetry point group P>2h and its eigenstates are con- 
veniently classified using the eight irreducible represen- 
tations of the dihedral point group. The bonding and 
anti-bonding eigenstates relative to the formation of the 
[2+2] cycloaddition belong to the B 3u and B lu symmetry, 
respectively (see Figure [5]) . 
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FIG. 5: Khon-Sham eigenvalues for the neutral system around 
the Fermi level (indicated by the horizontal line at e KS — 
e Fermi = ^)- The states are indicated by their corresponding 
symmetry labels (see text) and their energies as shown as 
a function of the intermolecular distance. Pictures at the 
inset (color online) show the density of the bonding (B3 U ) and 
anti-bonding (Bin) wavefunctions that undergo an occupation 
inversion at the transition distance (i?=8.90 A) 

The eigenfunctions of the Cgo are strongly modified 
when the intermolecular distance is changed, preventing 
us from studying the formation of covalent intermolecu- 
lar bonds as a function of the distance by a simple anal- 
ysis of the wavefunction close to the Fermi energy. The 
analysis of the Mulliken Overlap Population reported in 
the following section provides a useful tool for this pur- 
pose. The analysis of the eigenfunctions shows the exis- 
tence of four different regions of intermolecular distance 
where the eigenfunctions of a given symmetry exhibit the 
same ordering in terms of energy. The result is some- 
what different from that previously reported for a Cgo 



dimer by Fagerstrom and Stafstrom [nj. Figure [5] shows 
the HOMO-1, HOMO, LUMO and LUMO+1 eigenvalues 
and eigenstate symmetry for i?=8.5 A, 8.90 A, 9.35 A and 
10.85 A corresponding to R < R po i, R ~ R po i, R > R po i 
and R > R po i. The symmetry of the HOMO and LUMO 
eigenstates for increasing inter-center separation are A u , 
P>2u, Bi u and B211, B\ u , B 3u respectively. The presence 
of extra electrons results in a little modification of the 
general picture reported in Figure [5] mostly affecting the 
eigen energies. This suggests that for R < R po i the elec- 
tron doping would not affect the already formed cova- 
lent bonds since the B^ u and A g states at higher energy 
are non-bonding states. At R ~ R po i instead, the elec- 
tron doping would promote the occupation of the B\ u 
anti-bonding state, making of this a transition, unstable, 
state. At R > R po i and R ^ R po i additional electrons 
would fill the Bz u bonding state, contributing to favor 
the covalent bonds between the buckyballs. The pho- 
toexcitation of electrons from the anti-bonding B\ u state 
the bonding B 3u state would also favor the formation of 
covalent intermolecular bonds. The filling of this state 
is the reason for the decrease of optimum intermolecu- 
lar distance as a function of the extra charge, when the 
molecules are not yet polymerized. It is important to 
note that at these distances we observe a decrease of the 
HOMO-LUMO gap with the increase of R as shown in 
Figure O 



C. Charge Population Analysis 

The decomposition of the total charge into atomic con- 
tributions in molecules is not a uniquely defined prop- 
erty. Nonetheless, bond orders, or non diagonal elements 
of the density times the overlap matrix product provide 
useful information about the formation of intermolecular 
bonds. Here we choose to study Mulliken Overlap Popu- 
lations (MOP). The Mulliken analysis, although strongly 
dependent on the employed basis set (mostly through 
the overlap matrix in the product), is a reliable tool 
in cases where charge considerations are limited to the 
same atomic species. Figure [5] shows the evolution with 
the intersystcm distance and charge of the MOP between 
atoms 1-4 and 2-5 (see Figure[T|). Similarly to the analysis 
of the wavefunction symmetry, we have identified four re- 
gions with a different behavior of the MOP. At R < R po i , 
the MOP is around 80% of the value observed for the C-C 
single bonds in the Cgo molecule, consistent with the for- 
mation of covalent intermolecular bonds. The addition 
of electrons to the system does not modify this value 
because, as discussed in section IIII B[ these occupy the 
eigenstates with B^ u and A g symmetry corresponding to 
non-bonding states. When the center-to-center distance 
is close to the maximum of the energy barrier (8.88, 8.90 
and 8.95 A), the MOP discontinuously decreases, indicat- 
ing a weakening of the intermolecular bonds. We observe 
a value of ^0.05 electrons for the neutral system while 
smaller values arc found for the charged system. This 
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is due to the occupation of the anti-bonding B\ u eigen- 
state. At bigger intermolecular distances, at R > R po i, 
the polymerization does not occur since the MOP value 
drops to negative values indicating a repulsive overlap. 
The progressive occupation of the B^ u bonding state by 
1 and 2 electrons progressively increases the bond or- 
der of the cycloadduct and a positive overlap is observed 
(inset of Figure [5]) . The MOP value remains constant 
upon occupation of the B2„ state which is a non-bonding 
state whithin the four membered ring but it has a bond- 
ing character between atoms 3 and 6 and the other four 
equivalent pairs of atoms (see Figure [1]). This is the rea- 
son for the almost linear decrease of the second equi- 
librium intermolecular minimum for the unpolymerized 
fullerenes shown in Figure [2] A progressive occupation 
of bonding states between the molecules with increasing 
number of extra electrons will favor shorter intermolecu- 
lar distances and promote the [2+2]-cycloaddition reac- 
tion at lower pressures. 
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FIG. 6: Top panel: evolution as a function of the intermolec- 
ular coordinate of the Mulliken Overlap Population between 
atoms 1 and 4 (also 2 and 5, see Figure [TJ participating in 
the cycloaddition. The inset panel shows the zoom of MOP 
evolution for R > R po i. Bottom panel: evolution of the initial 
double C-C bond length between atoms 1 and 2 (also 4 and 
5). 



D. Structural Analysis 

Both buckyballs contained in the supercell are subject 
to an identical structural change with distance and elec- 
tron doping compatible with the P>ih symmetry. Such 
a change mostly affects the ten atoms belonging to the 
2 adjacent hexagons facing the neighboring molecule as 
shown in FigurclHland previously reported |20j . The bond 
length of the single bonds involving the atoms participat- 



ing in the cycloaddition (between atoms 1-4 and 2-5 in 
Figure [T|) slightly changes as function of distance and 
charge (up to a maximum of 0.08 A). On the other side, 
the initial double C-C bonds (between atoms 1-2 and 4- 
5) drastically change with distance (up to a maximum of 
0.20 A) and slightly with charge as reported in Figure [5] 
The result of the structural analysis is in agreement with 
the wavefunction and MOP analysis. At R < R po i, we 
observe a lengthening of the double C-C bond compared 
to the isolated molecule, up to a 10% increase, consistent 
with the formation of a [2+2]-cycloadduct. At these short 
intermolecular distances, the C-C bond between atoms 1- 
2 is not affected by the addition of electrons. At R > R po i 
instead, this C-C bond length decreases abruptly as ob- 
served for the MOP. In particular, it recovers the same 
value as in the isolated molecule showing that at these 
distances no polymerization of the molecules takes place. 
The addition of 1-2 electrons has a minor effect on the 
C-C bond length consistent with the MOP and wave- 
function analysis while the presence of 3-4 electrons does 
not further change this value. Instead, the C-C bond 
length between atoms 3-6 is reduced in agreement with 
the wavefunction analysis. 



IV. CONCLUSIONS 

We have performed ab initio DFT calculations of Cgo 
molecular chains in order to understand interplay be- 
tween charge doping and intermolecular separation on 
the polymerization process of the chains. 

Wc have studied the potential energy surface, the sym- 
metry of the wavefunctions around the Fermi level, the 
MOP between atoms of the cycloadduct and the evolu- 
tion of initial C-C double bonding length. 

The analysis of the potential energy surface shows that 
the total energy profile as a function of the intermolecular 
distance has a double-well shape. The MOP and struc- 
tural analysis show that the first and second minimum 
of the total energy profile correspond to polymerized 
and unpolymerized molecules, respectively. The equilib- 
rium intermolecular distance of the polymerized Ceo is 
not significantly shifted with charge doping. The [2+2]- 
cycloaddition reaction is prevented by an energy barrier 
which is related to the activation energy. The addition 
of electrons to the system results in an almost linear de- 
crease of the energy barrier as a function of the charge 
consistent with the study reported in Ref. [TH This is a 
very important result because it proves that in charged 
fullerenes the polymerization can occur at lower temper- 
ature. This is also consistent with the computational 
study performed by Sheka pol ] and the experimental ev- 
idence of a spontaneous one-dimensional polymerization 
for the alkali-intercalated fullerenes MCeo (M=K, Rb and 
Cs) under reversible solid-state transformation from the 
high-temperature phase 

The equilibrium intermolecular distance of the unpoly- 
merized molecules is appreciably affected by the electron 
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doping. In particular, the unpolymerized equilibrium dis- 
tance linearly decreases due to the progressive occupation 
of bonding wavef unctions. 

We conclude that polymerization depends on both the 
center-to-center distance of fullerenes and the negative 
doping of the system. Experimentally, pressure and tem- 
perature are required to approach the molecules and over- 
come the energy barrier, respectively. A different electron 
doping would affect the height of the energy barrier but 
not the distance at which polymerization occurs. In par- 
ticular, we observed that up to a doping of 4 electrons per 
2 molecules, the energy barrier progressively decreases. 
The size of the alkali cation clearly affects the intermolec- 
ular distance but also the charge transfer to the Ceo due 
to the decrease of the ionization potential with increas- 
ing ionic size. The potential energy curve suggests that if 
the intermolecular distance is R > R po i , for a given num- 



ber of metal atoms, bigger size cations should favor the 
polymerization due to the larger decrease of the energy 
barrier observed with higher electron doping. 



V. ACKNOWLEDGMENTS 

We acknowledge discussions with Xavier Blase, Al- 
berto Garcia, Enric Canadell and Berend Smit. MVFS 
was funded by DOE award number Ide-sc0003871 The 
work at the University of California were supported as 
part of the Center for Gas Separations Relevant to Clean 
Energy Technologies, an Energy Frontier Research Cen- 
ter funded by the U.S. Department of Energy, Office of 
Science, Office of Basic Energy Sciences under Award 
Number Ide-sc0001015l 



[1] A. M. Rao, P. Zhou, K. A. Wang, G. T. Hager, J. M. 
Holden, Y. Wang, W. T. Lee, X. X. Bi, P. C. Eklund, 
D. S. Cornett, et al., Science 259, 955 (1993). 

[2] Y. Iwasa, T. Arima, R. M. Fleming, T. Siegrist, O. Zhou, 
R. C. Haddon, L. J. Rothberg, K. B. Lyons, H. L. C. Jr., 
A. F. Hebard, et al., Science 264, 1570 (1994). 

[3] M. Nunez-Regueiro, L. Marques, J.-L. Hodeau, 
O. Bethoux, and M. Perroux, Phys. Rev. Lett. 74, 
278 (1995). 

[4] B. Sundqvist, Adv. in Phys. 48, 1 (1999). 

[5] I. B. Zhao, D. N. Poirier, R. J. Pechman, and J. H. 

Weaver, Appl. Phys. Lett. 64, 577 (1994). 
[6] S. Pekker, A. Janossy, L. Mihali, O. Chauvet, and 

L. Forro, Science 265, 1077 (1994). 
[7] P. W. Stephens, G. B. G. Faigel, M. Tegze, A. Janossy, 

S. Pekker, G. Oszlanyi, and L. Forro, Nature 370, 636 

(1994). 

[8] P. Launois, R. Moret, J. Hone, and A. Zettl, Physical 

Review Letters 81, 4420 (1998). 
[9] B. Verberck, H. Michel, and A. V. Nikolaev, J. Chem 

Phys. 116, 10462 (2002). 
[10] M. S. Dresselhaus, G. Dresselhaus, and P. C. Eklund, 

Science of Fullerenes and Carbon Nanotubes (Academic 

Press Inc.) (1996). 
[11] S. Margadonna, D. Pontiroli, M. Belli, M. R. T. Shikora, 

and M. Brunelli, J. A. Chem. Soc. 126, 15032 (2004). 
[12] E. A. Halevi, J. Phys. Org. Chem. 15, 519 (2002). 
[13] R. Poloni, D. Machon, M. V. Fernandez-Serra, S. Le 

Floch, S. Pascarelli, G. Montagnac, H. Cardon, and 

A. San- Miguel, Phys. Rev. B 77, 125413 (2008). 
[14] R. Poloni, G. Aquilanti, P. Toulemonde, S. Pascarelli, 

S. Le Floch, D. Machon, D. Martinez-Bianco, G. Morard, 

and A. San-Miguel, Phys. Rev. B 77, 205433 (2008). 
[15] R. Poloni, P. Toulemonde, D. Machon, S. Le Floch, 

S. Pascarelli, and A. San-Miguel, High Press. Res. 29, 

108 (2009). 

[16] M. R. Pederson and A. A. Quong, Phys. Rev. B 46, 13584 
(1992). 

[17] J. D. Santos, E. Longo, and C. A. Taft, Theo. Chem. 
625, 189 (2003). 



[18] J. Kiirti and K. Nemeth, Chem Phys. Lett. 256, 119 
(1996). 

[19] J. Fagerstrom and S. Stafstrom, Phys. Rev. B 53, 13150 
(1996). 

[20] E. F. Sheka, Chem. Phys. Lett. 438, 119 (2007). 

[21] A. S. Miguel, Chem. Soc. Rev. 35, 876 (2006). 

[22] J. M. Soler, E. Artacho, J. D. Gale, A. Garcia, J. Jun- 

quera, P. Ordejon, and D. S. Sanchez-Portal, J. Phys: 

Condens. Matter 14, 2745 (2002). 
[23] J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 

(1981). 

[24] C. Girard, P. Lambin, A. Dereux, and A. A. Lucas, Phys. 

Rev. B 49, 11425 (1994). 
[25] Y. Andersson and H. Rydberg, Physica Scripta 60, 211 

(1999). 

[26] J. Junquera, O. Paz, D. Sanchez-Portal, and E. Artacho, 

Phys. Rev. B 64, 235111 (2001). 
[27] E. Anglada, J. M. Soler, J. Junquera, and E. Artacho, 

Phys. Rev. B 66, 205101 (2002). 
[28] M. Leslie and M. J. Gillan, J. Phys. C 18, 973 (1985). 
[29] A. De Vita, M. J. Gillan, J. S. Lin, M. C. Payne, I. Stich, 

and L. J. Clarke, Phys. Rev. B 46, 12964 (1992). 
[30] R. Moret, P. Launois, P. A. Persson, and B. Sundqvist, 

Europhys. Lett. 40, 55 (1997). 
[31] G. Makov and M. C. Payne, Phys. Rev. B 51, 4014 

(1995). 

[32] S. H. Yang, C. L. Pettiette, J. Conceicao, O. Chesh- 
novsky, and R. E. Smalley, Chem. Phys. Lett. 139, 223 
(1987). 

[33] L. A. Girifalco, J. Phys. Chem. 96, 858 (1992). 

[34] C. S. Yannoni, R. D. Johnson, G. Meijer, D. S. Bethune, 

and J. R. Salem, J. Phys. Chem. 95, 9 (1991). 
[35] Y. Wang, J. M. Holden, X. X. Bi, and P. C. Eklund, 

Chem. Phys. Lett. 217, 413 (1994). 
[36] V. A. Davydov, L. S. Kashevarova, A. V. Rakhmanina, 

V. M. Senyavin, O. P. Pronina, N. N. Oleynikoy, V. Aga- 

fonov, R. Ceolin, H. Allouchi, and H. Szwarc, Chem. 

Phys. Lett. 333, 224 (2001). 



